Excitability in the vicinity of a saddle-node bifurcation: a mechanism for reversals. 
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We present a new mechanism for oscillatory or random reversals of the magnetic field that occur 
from a stationary dynamo state. The basic ingredient is the existence of two nearly critical dynamo 
eigenmodes, a stable and an unstable one, that collide during a saddle-node bifurcation. Above this 
bifurcation, a finite amplitude limit cycle is generated. In the neighbourhood of the bifurcation, even 
a small amount of turbulent fluctuations is enough to generate random reversals of the magnetic 
field. When the two competing eigenmodes have different symmetries, this scenario requires a broken 
symmetry of the flow that generates the dynamo. This is in good agreement with the observations of 
. the VKS experiment. It can also explain the dynamics of planetary or stellar magnetic fields and can 

be used to understand reversals of large scale flows that often develop on a turbulent background. 
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(N , INTRODUCTION 

f4 ■ Palaeomagnetic measurements show that the Earth magnetic field has a fixed polarity for long durations, but from 
Or time to time, it flips with the poles reversing sign (see the reviews by Dormy et al. (2000), Roberts & Galtzmaier 
Q ■ (2000)). The first simple models of field reversals involved disk dynamos (Rikitake 1958) or drastic truncations of 
D I the magnetohydrodynamic (MHD) modes (Nozieres 1978). A different class of models, also involving a few coupled 
differential equations, is based on the assumption that several magnetic eigenmodes are competing above the dynamo 
I threshold (Tobias et al. 1995, Knobloch & Landsberg 1996). In this framework, it has been proposed to relate 
. ' reversals to trajectories close to heteroclinic cycles that connect unstable fixed points ±B (Armbruster et al. 2001, 
Melbourne et al. 2001). Other models take into account turbulence in a stronger way by modelling it through random 
(— 1 ■ fluctuations of one of the coefficients of the dynamical system (Hoyng & Duistermaat 2004). A different approach, 
O ' initiated by Parker (1969), consists in trying to identify properties of the velocity field that generate dynamics of 
l— ~~ '■ the magnetic field. Recent numerical simulations have modelled the effect of the velocity with a fluctuating a-effect 
I (Giesecke et al. 2005, Stefani & Gerbeth 2005). Since 1995, three dimensional numerical simulations of the MHD 
■ equations in a rotating sphere have been able to simulate a magnetic field that displays reversals (see the reviews by 
£■ — ' Dormy et al. 2000, Roberts & Galtzmaier 2000, or for more recent works, Li et al. 2002, Kutzner & Christensen 2002, 
C*~) . Wicht & Olson 2004) . The experimental observation of reversals of the magnetic field and of other time dependent 
regimes has been performed only recently in a turbulent swirling flow of liquid sodium (VKS experiment) (Berhanu 
et al. 2007, Ravelet et al. 2008). 

The purpose of this paper is to give a new low dimensional description of the evolution of the magnetic field generated 
by dynamo action. From this description, different dynamical regimes of the magnetic field can be identified and their 
properties can be predicted. These predictions are in good agreement with the VKS experiment. They can be also 
tested against the results observed in numerical, geophysical and astrophysical dynamos. 



MODES AND SYMMETRIES 
Description of the VKS experiment. 



The VKS experiment involves a turbulent swirling flow of liquid sodium, generated by two impellers, counter- 
rotating at frequency -Fi (respectively F 2 ) in an inner copper cylinder, as sketched in fig. [1] (see Monchaux et al. 2007 
for the experimental set-up). When the disks counter-rotate with the same frequency F, a statistically stationary 
magnetic field is generated when F is large enough. Its mean value involves a dominant poloidal dipolar component, 
Bp, along the axis of rotation, together with a related azimuthal component B#, as displayed in fig. [1] (left). This 
shows that the VKS dynamo is not generated by the mean flow alone that would give a non axisymmetric magnetic 
field according to Cowling theorem. As explained in Petrelis et al. (2007), a possible mechanism is of a-cu type, the 
a-effect being related to the helical motion of the radially expelled fluid between two successive blades of the impellers, 
and the w-effect resulting from differential rotation due to counter-rotation of the impellers. 
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FIG. 1: Possible eigenmodes of the VKS experiment. The two disks counter-rotate with frequency F\ and F2. Left: magnetic 
dipolar mode. Right: magnetic quadrupolar mode. Poloidal (green) and toroidal (red) components are sketched. 

A striking observation of the VKS experiment is that time dependent magnetic fields are generated only when the 
impellers rotate at different speeds (Berhanu et al. 2007). When the rotation frequency difference / = F\ — F% is 
increased, the dynamo first remains stationary but the relative amplitudes of the field components are modified. The 
field is measured at position PI in fig. [T] A dipolar (respectively quadrupolar) field is associated to a field along B x 
(respectively B r ). Increasing the velocity difference from zero results in a continuous modification of the magnetic 
field from dipole to quadrupole. For a larger value of /, a bifurcation takes place and the field becomes time-periodic. 
The range of existence of this regime is small and if / is further increased, a stationary magnetic field is recovered. 
Remarkably, the stationary solutions are associated to points in the phase space that are located on the limit cycle 
related to the oscillatory solution. These bifurcations have been described in detail (Ravelet et al. 2008). It has been 
also shown that the parameter space (Fi,F 2 ) consists of several regions with an alternation of time dependent and 
stationary dynamos. Finally, a direct bifurcation from a non-dynamo state to a time dependent dynamo has been 
never observed. 

An amplitude equation for dipolar and quadrupolar modes 

As said above, the most striking feature of the VKS experiment is that time dependent magnetic fields are generated 
only when the impellers rotate at different speeds. We will show that this can be related to an additional symmetry. 
When F\ — F2, the experiment is invariant with respect to a rotation of -k around any axis located in the midplane 
between the two disks. Let TZ„ be one of these rotations. We expect that in the counter-rotating regime, the modes 
involved in the dynamics are either symmetric or antisymmetric. Such modes are displayed in figure [T] A dipole is 
changed to its opposite by TZ n , whereas a quadrupole is unchanged. More generally, we name "dipole" (respectively 
"quadrupole"), modes with dipolar (respectively quadrupolar) symmetry even though they might involve a more 
complex spatial structure. 

We assume that the magnetic field is the sum of a dipolar component with an amplitude D and a quadrupolar 
one, Q. We define A = D + iQ and we assume that an expansion in power of A and its complex conjugate A 
is pertinent close to threshold in order to obtain an evolution equation for both modes. Taking into account the 
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invariance B — > — B, i. e. A — > —A, we obtain 

i = fiA + vA + (3iA 3 + (3 2 A 2 A + (3 3 AA 2 + (3 4 A 3 , (1) 

where we limit the expansion to the lowest order nonlinearities. In the general case, the coefficients are complex and 
depend on the experimental parameters. 

Symmetry of the experiment with respect to TZ V when the disks exactly counter-rotate, amounts to constraints 
on the coefficients. Applying this transformation to the magnetic modes changes D into — D and Q into Q, thus 
A — > —A. We conclude that, in the case of exact counter-rotation, all the coefficients are real. When the frequency 
difference / is increased from zero, assuming that the coefficients are analytical functions of /, we obtain that the 
real parts of the coefficient are even and the imaginary parts are odd functions of /. 

When the coefficients are real, the growth rate of the dipolar component is fj, r + v r and that of the quadrupolar 
component is ji r — v r . The dipole being obtained for exact counterrotation implies that v r > for / = 0. By increasing 
/, we expect that v r changes sign and favors the quadrupolar mode. 

Study of the amplitude equation 

To analyse the properties of eq. [TJ we write the equations for the phase and the modulus r of A, A — r exp (i 0), 

r = r(n r + (3 2r r 2 + cos (20) (u r + {p Xr + (3 3r ) r 2 ) + cos (40) f3 4r r 2 

+ sin (20){vi + (#k - 0u) r 2 ) + sin (40) f3 4l r 2 ) , (2) 

= Mi + p 2l r 2 + sin (20) {-v r + {j3 lr - /3 3r ) r 2 ) - sm(40) (3 Ar r 2 

+ cos (20) ( Vi + (Jhi + fin) r 2 ) + cos (40) (3 4l r 2 . (3) 

The invariance B — > — B amounts to the invariance — > + tt. When the equation for does not have stationary 
solutions, the solution is oscillatory. Qualitative understanding is gained if we assume that the dynamics of r can be 
adiabatically eliminated. In this limit, r ~ 0, we obtain from eq. r ~ r o (0) such that we can write eq. [3] in the 
form 9 = Gi(0) - G 2 (0) and define 

Gi(9) = fi i + p2 i r 2 + cos(29)(v i + (p u + p 3i )r 2 ), 

G 2 (0) - sin (20) (v r -{(3 lr -/3 3r )r 2 ) + stn(49)f3 4r r 2 - cos (40) p 4l r 2 , (4) 

A fixed point C corresponds to an intersection of the curves G\ and G 2 - It is stable if G[(9 C ) — G' 2 (9 C ) is negative, 
i.e. Gi is larger than G 2 for smaller than C . It is unstable otherwise. 

When the solution is stationary, we obtain the relative importance of dipolar and quadrupolar components from 
D = ro cos(0 c ) and Q = sin(0 c ). For positive /j, r and v r , the equation has in general two pairs of solutions, one 
of which is stable and the other is unstable. It may have more solutions when non linearities associated to A 3 are 
important (Petrelis & Fauve 2008). 

When the values of the parameters are changed, the positions of the stable and unstable fixed points evolve. If they 
collide, they disappear and the solution becomes oscillatory. As an example, we consider the case = 1, (3 2r = — 1 , 
/ij = 1 and all the other coefficients are zero but v r . The functions G\ and G 2 are sketched in fig. [2] For fa < v r , 
there are 2x2 fixed points. If /ij > i> r , there are no fixed points and the solution is oscillatory. 

As said above, for the VKS experiment, we expect that u r changes sign if the frequency difference / increases from 
zero. When v r vanishes, if /i, is larger than z/,-, the solution is oscillatory. When v r reaches larger negative values, 
a stationary solution can become stable again. Therefore, we have a mechanism that leads to successive stationary, 
oscillatory and again stationary states when f is increased. 

Assume that the terms that vary, /ij, Vi and v r , are small compared to [i r which controls the amplitude of the fields. 
Then tq(9) weakly depends on 0. In the phase space (D = Re(A),Q — Im(A)), the stationary solutions are points 
located on the limit cycle associated to the oscillatory solution. 

To test this scenario, we numerically calculate the solutions of eq. Q]with the following coefficients: fj, r — 1; j3 2 = — 1; 
f3 1 = [3 3 = (3 4 = 0, v r = 0.05 (1 - 1.29 f 2 + 0.29 / 4 ), m = 0.0028 (/ + 2 / 3 ), Vi = 0.0104 (/ - 0.222 / 3 ), and / is a 
parameter that varies between and 1.5. Time series of the solution for various / are presented in fig. [2] (right). The 
phase space (Re(A), Im(A)) is shown in fig. [3] 
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FIG. 2: Left: functions G%(9) that control the phase dynamics defined by eq. [4] We have chosen fi r — 1. G\ (blue line), 
Gi (red line). The thin line corresponds to v r — 1.2. Stable fixed points (■) and unstable fixed points (•). The thick line 
corresponds to v r — 0.9, there is no longer any fixed point and the solution is oscillatory. Right: time series of the solutions of 
eq. Q]for / = 1.05 (thick line), / = 1.2 (thin line); dipolar component Re(A) (continuous blue line), quadrupolar component 
Im(A) (dashed red line). 
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FIG. 3: Phase space representation of the solutions of eq. [T] The stationary solutions are represented by symbols located on 
the limit cycle (obtained for / = 1) and are associated to increasing 6 when / increases from to 1.5 (/ = 0, 0.5, 0.75, 0.9, 
0.95, 1.2, 1.4, 1.475). 



When both modes are at the same distance from onset, v r ~ 0, the solution is oscillatory if /xf > vf + vf.. Therefore, 
when v r changes sign, the system evolves from a stationary solution towards an oscillatory solution and then back to 
a stationary solution. The limit cycle is generated by a saddle-node bifurcation and, slightly above onset, the slow 
phases of the cycle are very close to the former fixed points. This explains the form of the signal that displays a 
plateau or a slow evolution close to the former fixed points (see fig. [2] right). 

Appropriate scaling of the terms fj,i, Vi and v r as functions of / is required so that, the system is close to the 
quadrupolar solution before the bifurcation to the limit cycle takes place, the stationary solutions are points located 
on the limit cycle associated to the oscillatory solution, and when a stationary solution is recovered after the oscillatory 
regime, its phase 9 C increases for a large range of /. When these requirements are fulfilled, our model exactly reproduces 
the behaviour of the first bifurcations displayed by the VKS experiment when F\ — F 2 is increased from zero (see 




5 



figure 3 in Ravelet et al. (2008). 

SCENARIO OF THE BIFURCATION AND MECHANISM FOR REVERSALS 

A generic bifurcation 

In the VKS experiment, most dynamical regimes are reached through bifurcations from a stationary solution B s . 
At threshold, the generated limit cycle connects B s to — B s . It can be time-periodic, as it is the case for the regime 
described above, or random, and then looks like reversals of the Earth magnetic field. 

This kind of bifurcation is not classical (Arnold 1982). Actually, it is enforced by the symmetries of the problem 
and we can prove the following result. Consider a planar system invariant under the transformation B — > — B and 
with two different and non zero stationary solutions. One of the fixed points is unstable, B u and the other one is 
stable, B;,. The collision between the two fixed points generates a cycle that connects the collision point with its 
opposite, see fig. [4] This result can be understood as follows: the solution B = is unstable with respect to the 
two different fixed points, and their opposite. It is an unstable point, whereas one of the two bifurcating solutions 
is a stable point, a node, and the other is a saddle. If the saddle and the node collide, say at B c , what happens to 
initial conditions located close to these points? They cannot be attracted by B = which is unstable and they cannot 
reach other fixed points since they just disappeared. Therefore the trajectories describe a cycle. The associated orbit 
contains B = since, for a planar problem, in any orbit, there is a fixed point. Suppose that the orbit created from 
B c is different from the one created by —B c . These orbits being images by the transformation B — * B, they must 
intersect at some point. Of course, this is not possible for a planar system because it would violate the unicity of the 
solutions. Therefore, there is only one cycle that connects points close to B c and — B c . 

Regimes in the neighborhood of this bifurcation 

Two regimes can result from that scenario. One of the regimes is located above the onset of the saddle-node 
bifurcation and we now discuss its properties in the absence of fluctuations. The period T of the limit cycle diverges 
when the distance to the bifurcation threshold, e, vanishes. Indeed, for e very small, the period T is controlled by 
the time spent close to the former fixed points. There, the dynamics can be written, up to a change of variables, 
9 = 9 C + 4> (0 « 1), 

<j> = e + Ccj) 2 , (5) 

where C is a constant. This equation is the normal form of a saddle-node bifurcation (Arnold 1982). Integration of 
eq. [5] from the initial condition cf>(t = 0) = gives <f> = \f e/C tan (\feCt) from which the relation T ~ 2ir/\/7C is 
easily derived. 

The second regime is located below the bifurcation threshold and requires the presence of fluctuations. The system 
is attracted by one of the fixed points that is locally stable, ±B S . However, B s is very close to an unstable fixed point, 
B u (see fig. S]). Fluctuations can push the system to the unstable fixed point and once that point is reached, the 
system is attracted by the stable fixed point of opposite sign — B s . It remains close to — B s until fluctuations trigger 
a new evolution toward B s . This is a mechanism of excitability, that results in reversal-like behavior when the system 
is invariant with respect to B — > — B. If we approximate again the dynamics by the evolution of the phase <p, we can 
write, 

= e + C<j> 2 + C(t) , (6) 

where £ is a random noise of intensity D that describes the effect of the fluctuations. This equation is valid close to 
the bifurcation threshold for (f> <C 1 and we consider here e < 0. It can be written as the evolution equation of an 
overdamped motion in a metastable potential well V(<fi) = —e<fi — </> 3 /3. Between reversals, the system fluctuates in 
the vicinity of the stable fixed point. The durations of these "fixed polarity phases" are random. Classical results 
on the exit time from a metastable state show that they are exponentially distributed with a characteristic time 
that depends on the fluctuation intensity and the distance to bifurcation. Thus, P(T) oc exp (— Tj < T >), with 
< T >cx exp(AV/D) when the noise intensity D is small compared to the difference of energy between the fixed 
points AV = V((f> max ) — V(4>min) — 4C _1 / 2 (— e) 3 / 2 /3. This result is valid close to the saddle-node threshold so that 




FIG. 4: Phase space for reversals and excursions. The stable fixed points B a are black squares. The unstable fixed points 
B u are red circles. Initially (phase a), fluctuations drive the system away from the stable point (fluctuations driven dynamics 
sketched in black). Then, two situations can take place: in the first one, the system returns to B a (phase a', deterministic 
dynamics sketched in blue). In the second one, the system evolves to — Bs (phases b and c, deterministic dynamics in blue) 
and the system undergoes a reversal with an overshoot for the component Bi during phase c. 



eq. [S]is pertinent but not too close so that the fluctuations do not drive too frequent escapes. Nevertheless, this is the 
appropriate regime if we want to describe reversals that occur on time scales very long compared to the characteristic 
time scales of the deterministic dynamics. 

A model 

We use now a simple model in order to describe the effect of turbulent fluctuations on the dynamics of the two 
magnetic modes governed by [1] We change notation, A = B\ + iB 2 to stress that the two modes do not need to be a 
dipole and a quadrupole in general. 

B\ = + U r )Bi + {Vi - fli)B 2 

+C n Bf + Cai B\ B 2 + C 31 B 1 B\ + C 41 B\ + &nCi(*)£i + br 2 ( 2 (t)B 2 , 
B 2 = (jjir - v r )B 2 + (vi + fi t )Bi 

+C 12 Bf + C 22 B\ B 2 + C 32 B x B\ + C 42 B\ + br 3 <a(t)Bi + br 4 ( 4 (t)B 2 . (7) 

The nonlinear coefficients Cij are derived from those of eq. [T] Turbulent fluctuations are modelled by the terms Q that 
are independent gaussian white noises (with the Stratanovich interpretation). We take fi r = 1, v r = 0.05, \n = 0.045, 
Vi — 0.025 and (3 2r = — 1. The other non-linear coefficients are set to zero. For simplicity, we set br 2 = br 3 = and 
take br\ = br 4 = 0.1. Time series of the solutions are displayed in fig. O 

We observe that a reversal consists of two phases. In the first phase, the system evolves from the stable point B 3 
to the unstable point B u . The deterministic part of the dynamics acts against this evolution and the fluctuations are 
the motor of the dynamics. That phase is thus slow. In the second phase, the system evolves from B u to — B S) the 
deterministic part of the dynamics drives the system and this phase is faster. Both phases can be identified in fig. [5] 
(right). 

The behaviour of the system close to B s depends on the local flow. Close to the saddle-node bifurcation, the position 
of B s and B u defines the slow direction of the dynamics. If a component of B u is smaller than the corresponding 
one of -B s , that component displays an overshoot at the end of a reversal. In the opposite case, that component 
will increase at the beginning of a reversal. For instance, in the phase space sketched in figure [H the component B\ 
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FIG. 5: Left: time series of the solution of eq. [7]for bri — br^ — 0.1. Right: zoom on the first component. The horizontal lines 
represent the stable solution (upper line) and the unstable solution (lower line) in the absence of noise. 



decreases at the end of a reversal and the signal displays an overshoot. The component Bi increases just before a 
reversal. 

For some fluctuations, the second phase does not connect B u to — B s but to B s . It is an aborted reversal or an 
excursion in the context of the Earth dynamo. Note that during the initial phase, a reversal and an excursion are 
identical. In the second phase, the approaches to the stationary phase differ because the trajectory that links B u and 
B s is different form the trajectory that links B u and — B s . In particular, if the reversals display an overshoot this will 
not be the case of the excursion (see fig. [5] (right) and the sketch in fig. [4j. 

CONCLUSION 

We have studied dynamical regimes that can arise when two axisymmetric magnetic eigenmodes are coupled. 
Symmetry considerations allow to identify properties of the magnetic modes and, in some cases, put constraints on the 
coupling between the modes. We have shown that when a discrete symmetry is broken by the flow that generates the 
magnetic field, the coupling between an odd and an even magnetic mode (with respect to the symmetry) can generate 
a bifurcation from a stationary state to a periodic state. This behaviour is generic when a saddle-node bifurcation 
occurs in a system that is invariant under B — + — B. Close to the the bifurcation threshold, fluctuations drive the 
system into a state of random reversals that connect a solution B s to its opposite — B s . This scenario provides a simple 
explanation for many features of the dynamics of the magnetic field observed in the VKS experiment: alternation 
of stationary and time dependent regimes when a control parameter is varied, continuous transition from random 
reversals to time periodic ones, characteristic shapes of the time recordings of reversals versus excursions. 

Although the discrete symmetry involved for the flow in the Earth core is different from the one of the VKS 
experiment, a similar analysis can be performed for the geodynamo (Petrelis et al. 2008). 

More generally, our scenario can be applied to purely hydrodynamic systems. Cellular flows driven by thermal 
convection (Krishnamurti & Howard 1981) or by volumic forces (Sommeria 1986) display a transition for which a 
large scale circulation is generated on a smaller scale turbulent background. This large scale flow can display random 
reversals, very similar to the ones observed for the magnetic field. A model analogue to the present one, can explain 
how this large scale field can reverse without the need of a very energetic turbulent fluctuation acting coherently in 
the whole flow volume. 
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